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ABSTRACT 


In this study the entrainment-shear performance limit which occur™ 
in axial groove heat pipes was investigated and explained- In the exist- 
ing heat pipe literature the entrainment heat flux limit is defined as 
the condition where the Weber number is greater than or equal to one. In 
this analysis, the critical value for the entrainment Ueber number is 
found to be 2 iff _< We <_ 3it. Perhaps more important to the heat pipe de- 
signer than the entrainment performance limit is the prediction of the 
performance degradation due to vapor-liquid shearing stress which is also 
described. Preliminary qualitative experiments were conducted to observe 
the shear stress wave formation phenomena. Further more elaborate ex- 
periments would provide valuable verification of the shear stress wave 
phenomena and of associated performance degradation and limits. The equa- 
tions presented in this analysis may be used to predict and minimize the 
vapor-liquid shear stress performance effects that occur in axial groove 
and puddle flow artery heat pipes . 
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NOMENCLATURE 


2 

Vapor flow area, m 

2 

Liquid flow area, m 

2 

Averaged liquid flow area, m 
Wave velocity, m/s 

Wave velocity in absence of current, m/s 

Hydraulic diameter, 4A/P 

Internal pipe flow diameter, m 

Characteristic dimension, m 
Frequency , cycle/s 

Friction factor for liquid flow in groove 
Friction factor for liquid flow in circular tube 
Friction factor for vapor flow 

2 

Total vapor flow rate per pipe flow area, kg/m s 

Specific latent heat enthalpy, J/kg 

2 

Wick permeability, m 

Heat pipe effective length, m 

Characteristic mixing length, m 

Half-width of zone of influence at interface, m 

Wave number, 2 tt/A 

Number of pipe grooves 

Liquid transport factor 

. - 2 

Capillary pumping pressure rise, N/m 

.. - ■■ i - 2 

Pressure drop in the liquid, N/m 

- ■ •• - 2 

Pressure drop in the liquid due to viscosity, N/m 

Pressure drop in the liquid induced by counterflow, N/m 

... 2 

Pressure drop in the vapor, N/m 

2 - - ■■ 

Vapor pressure, N/m 

,-... 2 

Liquid pressure, N/m 

Wetted perimeter of the fluid flow passage, m 


Heat transfer rate, watts 
Effective pumping radius, m 

Reynolds number 

Dimensionless shear stress at vapor liquid interface, (St )/(u„U ) 

w A A 

Width of land between grooves, m 
Relative vapor-liquid velocity, m/s 

Axial vapor velocity, m/s 

Liquid velocity, m/s 

Groove width, m 
Weber number, dimensionless 
Coordinate parallel to mean flow, m 
Coordinate perpendicular to mean flow, m 


Greek letter nomenclature 


INVESTIGATION OP PERFORMANCE LIMITS 
IN AXIAL GROOVE HEAT PIPES 

I. 0 Introduction 

A survey of the various types of heat pipes used by NASA, ERDA, and 
industrial companies indicates a recent preference for configurations of 
axial and circumferential groove wick designs . ^ Such groove wick heat 
pipes , have the advantages of high evaporation and condensation film coef- 
ficients, reliable operation, ease of manufacturing, and potentially low 
cost. 

However, in operating the groove wick heat pipe, a number of per- 
formance limits are encountered. Since the groove wick heat pipe is nor- 
mally not used with additional screen wick, there is a significant counter- 
flow shear interaction between the vapor and liquid. The resulting pressure 
drop associated with this counterflow vapor-liquid shear results in lower 
than maximum heat transfer rates and in some cases unstable operation due 
to slugging or other entrainment-shear phenomena of the liquid. When 
operated in a gravity environment, the axial groove heat pipe may have a 
puddle of excess liquid accumulate along the lower side of the pipe. 
Counterflow shear and entrainment of liquid from the unshielded surface 
of the puddle can occur at relatively modest heat transfer rates. 

Also, while operating in a 1-g environment, the axial groove heat 
pipe can experience desaturation of the upper grooves along the top of the 
pipe. Depending on the construction of the grooves, these upper grooves 
may experience an additional hydrodynamic stress, which will cause them 
to first dry out- the upper side of the evaporator. Heat pipe performance 
is thereby restricted by this capillary limit in the grooves. Since these 
limits are not fully understood for axial groove or circumferential groove 
heat pipes , the expected heat flux limits and evaporation and condensation 
film coefficients needed for design of various groove geometries are not 
well known at this time. 

Because of these factors, a research project was initiated at The 
University of New Mexico to investigate the vapor-liquid shear and entrain- 
ment performance limits of the axial groove heat pipe. This project was 
supported by the NASA Ames Research Center on Contract No. NSG-2064 over 
the period from January 20, 1975 through January 19, 1976. 
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2.0 Previous Research 


A modest amount of previous rese?rch has been done on the problem of 
counterflow vapor-liquid entrainment-shear limits. A theoretical study 
was reported in 1968 by B&hr, Burck, and Huf schmidt 2 to predict friction 
factors and pressure drops associated with counter flow vapor-liquid shear. 
Also in this paper,- they described experimental work with water-air inter- 
actions which had been initiated. A companion paper published in 1969 by 

3 

Huf schmidt, Burck, DiCola, and Hoffman gives a detailed summary of the 
theoretical analysis results, and a comparison with the experimental data. 

4 

The theory was originally developed by DiCola assuming laminar flow 
of liquid in a rectangular open top channel with a counter flow vapor shear 
at the surface. The friction factor increase due to the vapor shear over 
that for Poiseuille of liquid in a tube was given by the ratio 


LP \ LP, 


fi ■ 


< 0.5 


( 2 . 1 ) 


6 T 


where S is the nondimensional stress, defined as 


w 




, and 


f LP / 2(l+i})) 2 (1-0. 627 <t>) 


( 2 . 2 ) 


is the friction factor ratio for S = 0, zero vapor shear, (t =0). 

For a heat pipe with N rectangular grooves of cross-section 6w and 
Ajj, = 6wN the interface shear stress was given by 


S = — S ' W — (f Re ) 
D 2 P V £ V V 


(2.3) 


where the symbols are defined in the nomenclature list. 
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Florschuetz has taken these results and applied them to the perform- 
ance of a sodium heat pipe and a water heat pipe both having open axial 
grooves. For both heat pipes the effect of vapor shear was found to re- 
duce the heat transfer rate significantly at low vapor pressures. 


A related paper by Zimmermann also deals with the counterflow shear 
friction factor and pressure drop in heat pipes. 


Another paper dealing with the counterflow shearing limit in the 

7 , " ' 

0 - g environment is by Schlitt, et al. The maximum heat transport is 

determined from the relationship 


Ap = Ap fl + Ap 
i v 

The liquid pressure drop consists of two terms, 

Ap A“ Ap £, o + Ap JMr 


(2.4) 


(2.5) 


where Ap„ is the viscous loss associated with the flowing liquid and 
Z,L 

applies to all capillary channels. The second term, Ap. , is the pres- 

J* f V 

sure drop in the liquid that is induced by the shearing action of the vapor 
counterflow and applies to open channels. 

The individual pressure drop ratios were given as 



Ap Ji,o ® L eff 

r 


ip s s £ 


#.! 

** ... •• 

Ap £,v ® L eff 

r 


Ap s 



Ap Q L ' r I 

*V * erf p 


’ _ V (w/D ) 2 1 

JL . y_ . f Re 

24 (1+s/w) , v v I 


Ap„ 




KA. 


1_ 

24 


V KA. 
v a 


_2 
D A 
v v 


f Re 
v v 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


where 


f v “ 


64/Re 


0 . 316/Re. 


1/4 


laminar vapor flow 
turbulent vapor flow 


(2.9) 


D. 


and K — 


2 

h,l 


lit, *e t ) 


(2.10) 


For the groove geometry given in Figure 1 , 
4A 

D h,£, 26+w 


( 2 . 11 ) 


and 


f £ Re^ = f(w/5) = 15 


( 2 . 12 ) 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 


The Equations 2.4, 2.5, 2.6 2.7, 2.8, and 2.9 can be combined to 
solve for the 0-g transport capability. The averaged liquid groove area 
was used which accounts for an average liquid meniscus recession. 



The presence of the f^Re^ terms in Equations 2.7 and 2.8 makes it diffi- 
cult to solve for the heat transport rate explicitly, particularly for 
turbulent flow where the Reynolds number has an exponent as given in 
Equation 2.9. A simpler formulation for predicting the shearing heat 
flux limit is needed for an explicit solution. 


. ’ SPACED o.Mfi»l«crrt ** «' f i* 

Fig. 1. Cross section of axial grooved 
aluminum extrusion. 

Also, a brief theoretical analysis of two-phase flow interactions 

8 

is given in the book by G. B. Wallis. However, none of these studies 
includes a thorough investigation of the entrainment interaction where 
waves and liquid droplets are formed. 

At the present time , in heat pipe design it is common to treat pre- 
diction of the entrainment based on prior knowledge of a "characteristic 

9 10 

dimension" in the heat pipe. ' Entrainment is presumed to occur when 
the inertial vapor forces exceeded the liquid surface tension forces at 
the vapor-liquid interface. The ratio of these forces is expressed by 
the Weber number 



p U 2 d 

We = (2.14) 

where d = characteristic dimension such as the space between wires in a 
screen wide heat pipe. 
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Cotter stated this xn a slightly different way. He stated that the 
wavelength A of a small disturbance on the liquid surface will grow in 
amplitude exponentially with time if, 

p U 2 A 

we = > 1 ’ (2 ' 15) 

"The process continues until the amplitude transcends the validity 
of the linear theory, and presumably ends in stripping of liquid 
droplets from the crests of the undulations. The quantity on the 
left is the Weber number with characteristic dimension A/2 tt, which 
is a geometric property of the capillary structure in the interface 
region." (Ref. 9, p. 347) 

Note here that Cotter's Weber number may be rearranged, 

2 

p U A 

We = - -- -- > 2 tt (2. IS) 

1,0 

Later Kernme' pointed out that the "entrainment heat flux limit," 

Q /A i is obtained by setting the Weber number equal to one , and combi n- 

Gfttl V 

ing it with 

• • • 

m = p U A and Q = m h^ , 
v v v fg 



(2.17) 


where 

P ev = mean vapor density at the evaporator exit. 

Equation 2.17 has been used to predict entrainment limit heat fluxes, 
by assuming a value for the characteristic dimension d. However, some in- 
vestigators have correlated this entrainment limit with d equal to the 
space between wires while others use the wire diameter. Obviously, there 
is a need for better understanding the entrainment phenomenon and to clarify 
the use of the characteristic dimension d. 
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3.0 WAVE FORMATION AN!? ENTRAINMENT 


Since the previous research just described has defined the friction 
factors sui cable for use in predicting increased pressure drop in axial 
groove heat pipes but has not considered the unstable effects of waves, 
ripples, or entrainment droplet formation, the work, presented in this . 
report is focused primarily on the latter phenomena. Also, a discussion 
is presented to help clarify the use of the wave length, A, the character- 
istic dimension, d, the Weber number, and the entrainment-shear heat trans- 
fer rate limit. 

11,12 

Hanratty, et al. have conducted both theoretical and experimental 

studies on the initiation of waves due to the flow of gas or vapor over liquid . 
Parts of the results of these studies have been incorporated in the follow- 
ing discussion. Suppose a liquid is flowing in one direction in a gently 
downward sloping open channel and its vapor is flowing over it in the op- 
posite direction, as shown in Figure 2. At zero or very low vapor velocity, 
the flowing liquid will be undisturbed, however , with increasing relative 
velocity, a destabilizing effect is noted. The first waves to appear are 
small two-dimensional ripples which extend over the width of the channel 
and travel in the direction of the vapor flow. At even higher vapor velo- 
city the amplitude and the wave velocity increases and a three-dimensional 
cross-hatched wave pattern appears. This latter case is similar to the 
waves which are obtained on a body of water during a light squall. 

At a vapor velocity which is about double that necessary to produce 

the cross-hatched waves, the first roll waves appear. The roll waves 

occur when the crests of the waves are moving faster than the troughs . 

This causes the downstream end of a wave to steepen and to roll over upon 

itself. Thus, liquid in the wave is moving in a circular path. Such waves 

13 

are called roll waves. They were first described by Cornish who ob- 
served them in mountain streams. A good review of the literature on this 
. 12 

topic is given by Hanratty. The first theoretical prediction for the 
conditions causing roll waves was made by Jeffreys 14 who considered the 
case of turbulent flow in an open inclined channel. 

At higher vapor velocity the shear forces on the tops of the roll 
waves and the natural rolling action of the waves are adequate to over- 
come the surface tension forces of the liquid so droplets of liquid are 
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a) Undisturbed liquid b) Drag on liquid initiates 

flow reversal at surface, 
and ripple waves are formed. 




e) Initiation of droplet 
entrainment and "white 
water . " 


Fig. 2. Subsequent wave phenomena 
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REPRODUCIBILITY OF THE 
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separated from the surface and entrained into the vapor stream. The onset 
of droplet entrainment is usually preceded by a noticeable roughening of 
the tops of the roll waves which resembles patches of "white water." 
Further increases of vapor flow velocity increase the entrainment. 

The sequence of events in the entrainment process is shown in Figure 
2. The corresponding regions of wave activity are shown in Figure 3 for 
an experiment conducted by Wallis, et al. in which air flowed over water 
in a horizontal duct 30 cm (12 inches) wide and 12.7 cm (5 inches) deep."*" 5 

In a heat pipe the onset of entrainment may be encouraged by nucleate 
boiling in the evaporator which can occur at sufficiently high heat fluxes 
The bursting of bubbles at the liquid surface, as shown in Figure 4, has 
long been known as a source of droplet entrainment. 

In order to help elucidate the phenomena of wave formation and en- 
trainment, the following definitions of terms are givens 

1. Ripple Waves: Two-dimensional waves formed by gentle undula- 

tions of low amplitude (less than 0.1 mm, 0.005 inches) on the 
liquid surface which travel at about 22-30 cm/sec (0.75 to 1.00 
feet/sec) in an air-water system, and are not initiated by tur- 
bulence in the vapor flow. 

2. Cross-hatched Waves: As the vapor velocity is increased 

slightly these three-dimensional waves appear as a criss- 
cross ripple or as a pebbled pattern. For an air-water 
system the wave lengths are of the order of 0.5-1 cm (0.25- 
0.5 inches) and the waves travel at only slightly greater 
velocity than the ripple waves. 

3. Roll Waves: As the vapor velocity is increased the cross- 

hatched waves become stable and persist up to a vapor velo- 
city about double that necessary to produce the cross- 
hatched waves. At this point, roll waves appear and travel 
at about 60 cm/sec (2 ft/sec) down the channel. The roll 
waves may take a finite distance along the channel to form 
and normally they appear first where the film depth is a 
maximum. 
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Droplet Entrainment: Further increases of vapor velocity, 

up to about 50 percent more than that needed to produce roll 
waves, will cause droplets to be torn from the liquid sur- 
face and entrained into the vapor flow. 
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4.0 EXPERIMENTAL OBSERVATIONS 


In order to observe the entrainment-shear interactions, three dif- 
ferent experimental heat pipes were built s -.nd preliminary qualitative ob- 
servation tests were conducted. A schematic drawing of the final appara- 
tus is shown in Figure 5 and a photo is shown in Figure 6._ The apparatus 
is composed of an axial grooved copper tube with a 91.4 cm (3 foot) evapo- 
rator and condenser section, and a 15.24 cm (6 inch) glass tube adiabatic 
transport section. Each end of the tube has glass observation windows for 
viewing phenomena along the axis. 

In the experimental tests Freon-12 was used as the working fluid 
since it has a low latent heat of vaporization which will give large mass 
flow rates at modest heat transfer rates. With the heat pipe in a hori- 
zontal orientation, the evaporator would ejqperience nucleate boiling and 
dryout before any wave phenomena were observed in the adiabatic transport 
section. By raising the condenser end 1.27 cm (1/2 inch), and by adding 
cooling water over the condenser to lower the pipe temper .ure to 27 °C 
(80°F) , waves were observed in the adiabatic transport section, as shown 
in Figure 7. The average wavelength of the waves observed was approxi- 
mately 0.4 cm (0.157 inches). Also, in Figure 8, a sketch of the observed 
wave phenomena is shown. The waves gradually' increase in amplitude down 
the tube from the boiling wave source in the evaporator and become fully 
developed in the adiabatic transport section and gradually die down as 
they enter the condenser section. Further increasing of the condenser 
elevation, evaporator heat input rate, and a slight reduction of the con- 
denser temperature showed a slight -increase in wave speed but no entrain- 
ment was observed. Insufficient time and equipment prevented further 
qualitative experiments. A chilled liquid cooling unit was ordered, but 
was not received in time to conduct tests at lower temperatures . One 
would expect that some entrainment could be observed at lower temperatures 
where the vapor pressures and densities are lower and the entrainment limit 
is lower . 
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1.58 cm (5/8 in.) OD 
Axial Grooved Copper 
Tube 

Glass Tube in 
V / Midsection 


Condenser 


Electrically 

heated 


Glass window 
on each end 


Cooling 

Water 


Fig. 5. Sketch of experimental apparatus 


reproducibility of th 
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Photograph of experimental apparatus 





Fig. 7. Photograph of waves observed in the adiabatic 

transport section. (The evaporator is at the right 
hand side of the photo.) 



Fig. 8. Sketch of wave phenomena observed 




5.0 HYDRODYNAMIC WAVE INSTABILITY 

The hydrodynamic instability of waves formed on a liquid surface by 
gas or vapor flow over it may be explained physically as follows . As the 
surface film is disturbed from its equilibrium position, normal and tan- 
gential stresses are produced and waves will form. The normal stress com- 
ponent is in phase with the wave displacement, and it attempts to deform 
the liquid surface by exerting an upward force on the crests and a down- 
ward force on the troughs of small periodic waves. Also, the tangential 
stress component is in phase with the wave slope, and it has the effect 
of accelerating the liquid on the windward slope while retarding that on 
the leeward slope. The latter mechanism tends to displace fluid towards 
the crests and away from the troughs of the waves, thereby increasing 
their amplitude. 

The crests of the long wavelength disturbances on a flowing liquid 
move faster than their troughs and this causes the downstream part of the 
wave to steepen and roll over upon itself. The thinner the liquid film, 
the more effective this mechanism is, for in very thin films, the pertur- 
bation velocity in the horizontal direction is much greater than in the 
vertical, so the influence of the tangential and vertical stresses extends 
throughout the whole film. The development of waves is governed by con- 
tinuity conditions . The wave growth rate is proportional to the rate of 
accumulation of liquid near the crests, and this, in turn, is inversely 
proportional to the liquid viscosity. 

Wave instabilities grow when the disturbing vapQE - stress is sufficient 
to overcome the restoring forces of gravity and surface tension. With a 
given vapor flow rate, instabilities can always be induced in an initially 
stable film by decreasing its thickness. Also, increasing the liquid flow 
velocity simply increases the roll wave frequency. 


5.1 Analysis of the Wave Instability 

As is shown in figure 9 , the wave flow along the surface of a liquid 
film caused by a vapor flow may be considered to be a sinusoidal wave mov- 
ing at speed C. Analysis of the problem is simplified by using a coordinate 
system moving at the speed U^, so that the liquid flow velocity, U^, is zero 
and the relative vapor velocity is = U £+■' 


y 


Core (Negligible liquid) 



Pig. 9. System for film stability analysis. 
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The dynamic equation for the surface wave given by Tippets is 

% - e 2 n Q = o ■ (5.1) 

where n Q is the amplitude of the surface wave as shown in Figure 9 and 3 
is a disturbance growth parameter given by 


3 » 


(us-mO^) ^ tanh(m&' 1 ) + — — m 2 coth(m6) 

P v 


- am /p 


Pp 

— coth(m6) + tanh(m£," ) 


(5.2) 


and the radian frequency of the interfacial waves is 
mU 


U) = 


r£ 


P £ coth(m6) 
P v tanh(m£") 


(5.3) 


2 

The vapor- liquid interface will grow indefinitely with time if g > 0. 
The critical wavelength for 3 = 0 , is given by 



(5.4) 


The fastest growing disturber ;es are those having a wavelength X which 

m 

maximizes the growth factor 3> which is given by 
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x - 

m 


3 TTQ 

P u 2 g 
v vH 




(5.5) 


In obtaining these 


results Tippets 
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assumed that the fluid is in 


potential flow. This idealization approximates the condition of well de- 
veloped turbulence at the interface. He related the wavelength and film 
thickness at instability to the common properties of the flow, under the 


assumption that the interfacial wave motion is dominated by the stream 
turbulence. He used Prandtl's mixing length theory to relate the turbu- 
lent shear stress x = -pu'v' to the local mean flow properties by a char- 
acteristic mixing length This analysis assumes that the surface wave- 

length, X, is proportional to & m , V' , and to the film thickness, 5. He 
further explained that due to the character of turbulent flow there exists 
a spectrum of possible disturbance wavelengths extending from near the dimen 


sions of the molecular motion up to within an order of magnitude of the 
duct dimensions. Disturbances with wavelengths A < A^ can be ignored, 
since their amplitude will remain small. Disturbances with wavelengths 
X > A^ can also be ignored, because their growth rates will be no faster 
than those with wavelengths slightly less than A^. Thus, the dominant 


wavelength A q of the characteristic disturbance motion impressed on the 

interface by the turbulence is bounded according to A < A < A . Hence, 
J c o m 

from Equations 5 . 4 and 5.5 



( 5 . 6 .) 


The corresponding film thickness is 


6 


<V*i? 


P v U v£ 



(5.7) 


where is def ined by A q = K^<5 , and ~ 1.0. 

The entrainment begins after the initiation of the instability if the 
disturbance wavelength is within the A^ range. The occurrance of entrain- 
ment reduces the liquid film thickness'. When the film; thickness is re- 
duced the wave disturoances are reduced, thereby reducing the entrainment. 
Thus, the wavelength and film thickness will tend to stabilize about the 
values given by Equations 5.6 and 5.7. 
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Heat transfer by evaporation from the liquid layer may be included in 
the above analysis. However, boiling heat transfer cannot be included 
since the forces on the film associated with boiling bubbles have not been 
incorporated in the analysis. It would be expected that liquid films on 
high heat flux surfaces where boiling occurs would tend to be thinner than 
that predicted by Equation 5.7. 

5.1.1 Weber Number Relationship 

We can write the Weber number, previously given by Equation 2.14, as 


We 


p u 2 .a 
v vZ o 


(5.8) 


where X is used for the characteristic dimension, in place of d. We now 
o 

note that if Equation 5.6 is arranged as 


^-T^= K 3 (1 + Pv/PA> 


then 


We = K 3 ( 1 + p v /p £ ) , 2ir < K 3 < 3 tt. (5.9) 

If we can assume that p v «p^ (for example, with water at 100°C, p^/p^ = 
6.24 x 10 *V, then 


We = K 3 , 2 tt < K 3 < 3tt. (5.10) 

Thus, according to this analysis, growing instabilities and entrain- 
ment will begin when We > 2r. rather than We = 1 as has been suggested in 

9 

the past. However, if Cotter's Weber number is interpreted as in Equa- 
tion 2.16, then there is agreement that growing instabilities and entrain- 
ment occurs when We > 2 t. 

5.1.2 Wave Speed Relationship 

The initiation of wave instabilities may be further described by de- 
termining the wave speed and wave length at the onset of ripple waves. . The 
wave speed, C, may be most easily determined by a coordinate system with 
the x-axis at the averaged liquid -vapor interface. Again, a sinusoidal 
wave is assumed to exist at the liquid-vapor interface as shown in Figure 
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The analysis, given in Appendix A, describes the wave speed by the follow- 
ing relationships. 


Vi + p v p v> ± J _ Vt tu v + °t l 


p l • f <V 


(!> t + M 


1/2 


(5.11) 


where 


c 2 = A 


(p £ ~ p v ) 2m q 

2 tt (p A + P v ) X (P„ + Pj 


(5.12) 


£ v 


.. i s - 


C Q is bhe wave speed of a traveling ripple wave in a quiet fluid where 
U = U- = 0, that is, like that observed oy a pebble thrown in a still pond. 

V Xi ■■■■" - :: 

From Equation 5.11, we can see that the wave will be unstable if 

2 , 1 . . . ....... . .. , 

> C 2 (5.13) 


The value - of X for which C q is a minimum may be designated by X . 

and corresponds to the point of onset of ripple waves. This value of X - - 

.. ” " rip 

can be determined by differentiating Equation 5.11 with respect to A and 

setting it equal to zero. We then find that the value of X . is 

.... rip 


p £ p v 


u 


v£ 


P £ + P v 
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* 





and the corresponding minimum value of C q is 


C 

c 


rip 


20 (Pa - » v ) 

<p t + p v> 2 



(5.14) 


(5.15) 


For a flowing vapor or gas to cause waves on a liquid surface, its 

speed must be greater than or equal to the minimum ripple wave speed, that 

is, U . > C . As the vapor velocity is increased, the wave speed and 
v* °rip 

amplitude will increase correspondingly, decreasing in wavelength until 

A ~A (where K_ - 3ir) is reached, where we assumed in the analysis that 
o m 

the entrainment is well developed. 


5,1.3 Unstable Wavelength. 

From Equation 5,13, the liquid film will be unstable if 

Vv [ » c o • 


(5.13) 


and we know there is a dominant wavelength given by Equation 5. 6 


A = 
o 




2ir < < 3 it 


(5.6) 


v vX 


after rearranging and dividing through both sides by (p + p ) and dropping 

a v 

the A q subscript, we get 


P„X 


vX 


X v l p„ + p 


gK 3 

Mp 4 + p„> 


(5.16) 


The righthand side of Equation 4.16 gives the value of P^P V 
term of the wave length A. 


( \ 2 ■ 
VP^) in 


Thus, the incipient entrainment wavelength can be found by using Equa- 
tions 5.12 and 5.16 satisfying the relationship 5.13. This can be done by 

■ 2 
plotting C q vs . A and 
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vs. X 


3 

X(P A + p v ) 

the intersection point between two curves gives the value of incipient 

entrainment wavelength, X , as shown in figures 11 and 12 for ammonia 

en 

at 203 K and 320 K respectively. 
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Entrainment wavelength for ammonia, T = 203 



Entrainment wavelength for ammonia, T = 320 




6.0 ENTRAINMENT HEAT FLUX 

The entrainment is assumed to be developed as the result of the 
liquid surface dynamic instability at wavelength 


X 

en 


( i + M 

v ' * 1 


thus, 


D vl = 


3tto 


P X 
v en 



( 6 . 1 ) 


( 6 . 2 ) 


For relatively low liquid velocities the following approximation holds, 


U 


v*. 



(6.3) 


and because 



U 

v 



so 


SL 

A 


) ■ 
en 


P U . h. 
v vf fg 


(6.4) 


By combining Equations 6.2 and 6.4, the incipient entrainment heat flux 
can then be determined by, 



(6.5) 


Equation 6.5 has been plotted along with the sonic heat flux limit and the 

wicking limit heat flux including counterflow shear effects for water , 

freon-1.1, ammonia, and methane in Figures 13 through 16, respectively. 

18 

The equation used for the sonic heat flux limit is that given by Levy 


• \ P„ h- 

« y -5 £ g 

A v / 2 /”k+l~ 

f e 


( 6 . 6 ) 


where p and C s are the vapor density and speed of sound evaluated -at the 
evaporator entrance temperature and k = c /c is the ratio of specific 
heats. Equations 2.4 through 2.13 given by Schlitt, et al.^ were combined 
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Heat Flux, Q/A , watts/m 



Fig. 13. Performance heat flux limits for water in an axial groove 
heat pipe. 
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Temperature, K 


Fig. 14. Performance heat flux limits for freon-11 in an axial 
groove heat pipe 










Heat flux, Q/f\ , watts/m 



Fig. 15. Performance neat flix limits for arroonia in an axial 
groove heat pipe 
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Temperature, K 


Fig. 16. Performance heat flux limits for methane in an axial 
groove heat pipe 










to give, 


-1 


Q w 


r p L eff 


8 \ 


(w/D v ) 


3 V (1+s/w) 


+ 1 + 


8 V v 


KA„ 


3 d 2 a 
v v 


(6.7) 


which gives the wicking limit heat transfer rate including vapor- liquid 
counter flow shear effects . The axial groove heat pipe dimensions given 
in Figure 1 were used in these calculations . 

In the Figures 13 through 16 it is apparent that the sonic limit heat 
flux is considerably higher than the other limits for all the fluids. How- 
ever, since the wicking limit is more sensitive to the liquid properties 
than the other limits considered here, the wicking limit exceeds the en- 
trainment limit for water as shown in Figure 13. For freon-11, ammonia, 
and methane, which are less effective working fluids (N^ is lower) , the 
wicking limit is the lowest limit. In all cases, the limits are lowest 
at low temperature where vapor densities are low and vapor velocities are 
high. This low temperature condition occurs during startup and causes re- 
latively large pressure losses and low heat transfer rates. 
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7.0 CONCLUSIONS AND RECOMMENDATIONS 

According to the analysis presented in this report, using the lin- 
earized instability theory, wave instabilities grow and entrainment begins 

when 2 tt < We < 3 tt rather than at We > 1 as suggested in the past. Here 

2 

the Weber number is defined as We = (p U e X )/a and X , which is the wave- 

V vx> o o 

length of the growing instable wave defined by Equation 5.6, is used as 
the characteristic dimension in place of d. It should be noted that this 
result is identical to that suggested by Cotter in Equation 2.15 if both 
sides of his equation are multiplied by 2t\ to give Equation 2.16. Cotter 
further stated that X/2 tt = d, is a characteristic dimension which is a 
geometric property of the wick surface. Additional experimental data is 
needed in order to determine any relationship between d and the wick sur- 
face geometry, particularly for axial groove heat pipes. 

The effect of the counter flow shear force between the vapor and 
liquid is to create an additional pressure drop, given by Ap^ v in Equa- 
tion 2.7, which reduces the overall heat transfer rate of the heat pipe, 
given by Equation 6.7. If wicking forces are sufficiently high, as they 
are for water, the wicking limit can exceed the entrainment limit which 
was shown in figure 13. Thus, according to the calculations plotted in 
figure 13 for the axial groove heat pipe shown in figure 1, it is possi- 
ble to observe entrainment of water since the entrainment limit occurs . 
before, the wicking limit. However, for this axial groove heat pipe ope- 
rating with ammonia, freon-11, or methane, as shown in figures 14, 15, 
and 16, the wicking limit would be reached before the entrainment limit 
and entrainment would never be attained. 

The entrainment limit heat flux is given by Equation 6.5, the entrain- 
ment velocity by Equation 6.2, and the entrainment wavelength by Equation 
6.1. The onset of ripple waves occurs at a wave speed and wavelength given 
by Equations 5.15 and 5.14, respectively. 

Further experimental measurements of the counterflow shear effects 
and the entrainment limits need to be made and compared to the values pre- 
dicted by the equations given in this report. 


8.0 


References 


1. K. T. Feldman, Jr. , Editor, Heat Pipe Technology Bibliography , The 
University of New Mexico, Albuquerque, 1974. 

2. A. Bahr, E. Burck, and W. Huf schmidt, "Liquid-Vapor Interaction and 
Evaporation in Heat Pipes , " Paper D-6 , 2nd International Conference 
on Thermionic Electric Power Generation, May 27-31, 1968. 

3. W. Hufschmidt, E. Burck, G. DiCola, and H. Hoffman, "Der EinfluuS 
der Scherwirkung des Dampfstromes auf den laminaren Flilssigkeitsstrom 
in Kapillaran von WKrmer ohren , " Wiirme-und stof f ilbertragung , Bd. 2 
(1969) S. 222-239. 

4. G. DiCola, "Soluzione Analitica, a mezzo della transformata di Fourier 
di un problema flusso in a canale rettangolare , " EURATOM, C.C.R. , 

Ispra, Italy, C.E.T.I.S., 1968. 

5. L. W. Florschuetz, "Estimates of the Significance of Vapor-Liquid 
Shear Interaction on the Performance of Heat Pipes with Open Rec- 
tangular Groove Capillary Structures," Los Alamos Scientific Labora- 
tory, N-5 Technical Memorandum No. 172, Los Alamos, New Mexico, 

August 10, 1970. 

6. P. Zimmermann , "Einfluss der Dampfstr&mung auf das instationSre 
Betriebsverhalten von warmerohren , " Dissertation, Universit&t 
Stuttgard, (1972) 

• « 

7. K. R. Schlitt , J. P. Kirkpatrick, and p. J. Brennan, "Parametric Per- 
formance of Extruded Axial Grooved Heat Pipes from 100° to 300°K, " 

AIAA Paper No. 74-724, Thermophysics and Heat Transfer Conference, 
Boston, Massachusetts, July 15-17,,. 1974. 

8. G. B. Wallis, One-Dimensional Two-Phase Flow , McGraw-Hill Book Company, 
1969. 

9. T. P. Cotter, "Heat Pipe Startup Dynamics," 1967 IEEE Thermionic Con- 
version Specialist Conference, Palo Also, California, October 30 - 
November 1, 1967, pp. 344-348.- 

10. J. E. Kemme, "High Performance Heat Pipes," 1967 IEEE Thermionic Con- 
version Specialist Conference, Palo Alto, California, October 30 - 
November 1, 1967, pp. 355-358. 

11. T. J. Hanratty and J. M. Engen, "Interaction Between a Turbulent Air 

Stream and a Moving Water Surface," AIChE Journal, Vol. 3, September 
1957. " ■ “ - - . ■ : .. 

12. T. J. Hanratty and A. Hershamn, "Initiation of Roll Waves," AIChE - 
Journal, Vol. 7, September 1961. 

13. Vaughn Cornish, Waves of the Sea and Other Water Waves , T. Fisher 
Unwin, London, England, 1910. 


30 


14. H. Jeffreys, Proceedings of the Royal Society, A107, 189, 1925. 

15. G. B. Wallis/ J. Mr Turner, I. Beiriberis, and D. Kaufman, AEC Report 
NYO-3114-4, 1964. 

16. G. F. Hewitt, N. S. Hall-Taylor , Annular Two-Phase Flow , Fergamon 

Press , 1970. > 

17. F. E. Tippets, ^Analysis of the Critical Heat-Flux Condition in High- 
Pressure Boiling Water Flows," Journal of Heat Transfer, Trans. ASME, 
February 1964, pp. 23-38. 

18. E. K. Levy, "Theoretical Investigation of Heat Pipes Operating at Low 
Vapor Pressures," ASME Transactions - Journal of Engineering for 
Industry, November 1968, pp. 547-552. 

19. N. E. Kochin, I. A. Kible, N. V. Roze, Theoretical Hydrodynamics , 

John Wiley,. 1964. -- . . ' 


31 




APPENDIX A 


Detailed Analysis of Wave Instability 

Here an original analysis of the counterflow shear wave instability 
is presented. For this analysis a fixed coordinate system is assumed 
where a liquid film is moving at '(peed u^ and a vapor moving over it at 
speed u as shown in Figure Al. 



The liquid and vapor flows may be represented by a velocity potential, 

<j> = ux, plus a perturbation term, <f>' , to account for the waves, 

A = u x + 4>' (A. la) 

r v v 

<t> z = + <j>' z {A. lb) 

The corresponding boundary conditions are 

<j>' = 0 at y = + i" 
v 

and 

<p' z = 0 at y - - Z" 

The use of velocity potentials requires that the flow is irrot .tional,- 
that is, the perturbation terms, tj > ' , must satisfy the Laplace's eq;ation, 


(A. 2a) 
(A. 2b) 
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+ aft „ o 
2 2 
Bx By* 


(A. 3) 


The perturbation terms may then be described by a sinusoidal term, 

2 TT 


<f>' (x,y,t) = Y(y) sin ■— (x - Ct) 
Introducing A. 4 into A. 3 we get 


(A* 4) 


f§ 

ay 


where the solution is 


2 7T 


2ir 


A y “ t~ y 

Y » Ae + Be 

thus Equation a. 4 becomes 


2ir 


2ir 


y — — y 

^(XrY/t) - (Ae + Be * ) sin — (x-Ct) 


(A. 5) 


Applying the boundary condition A. 2a to ' we get 


4>Wx,y,t) - D sin -p- ( y-V ') sin ~ (x-ct) 


(A. 6a) 


and the boundary condition A. 2b gives 


(x,y,t) = E sinh (y+i") sin ~ (x-Ct) 


(A. 6b) 


where D and E are constants. 


The forms of Equations A. 6 indicate that the displacement rs has a 
sinusoidal form as. 


2tt 

r) - t] cos ~ (x-Ct) 
O A 


(A. 7) 


The linearized kinematic conditions across the liquid-vapor inter- 
face are 
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3ri 3H. _ __i 

9t U v 3x 3y 


(A. 8c;) 


3jX ♦an v ii 
at U £ 3x 3y 


(A. 8b) 


Substituting Equations A. 6a, A. 6b, and A. 7 into A. 8a and A. 8b respect- 
ively, we get 


n (C-u ) 
o v 

cosh <y-£") 


(A. 9a) 


n o (C ~V 

cosh 1^ ty+£") 


(A. 9b) 


The dynamic pressure gradient condition across the interface in 
linearized form is 


n - K - - “ fi 

3x 


(A. 10) 


where p 1 is the pressure perturbation and a is the surface tension. 

The pressure, p* , may be calculated from the linearized Bernoulli 
equation 


r “ >t + u aT ' 99 


CA.ll) 


Substituting Equation A. 11 into A. 10, we get 


‘ P V + U V aT ' gn) p » ( 3T *"137 ' 9 ’' ) 


(A. 12) 


Now we substitute Equations a. 6a, A. 6b and A. 7 into A. 12 and eliminate 
D using A. 9a, E using A. 9b, and n Q , we get 
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p (-C+u ) (C-u )tanh v— (y-£") 

V V V A 

p £ (-C+u £ > (C-u £ )tanh {y+£") - (P £ - p v ) = o ~ 

The above equation may be simplified with very little error, as did 
Tippits fl7] by recognizing that £" = .4 A q and y is very small so that 
the tanh terms are nearly equal to 1, as 


P v (C-u v ) 2 + P £ (C-u £ ) 2 - Si 


/ \ 2-rra _ 

2ir P £ P v ) “ A ~ ° 


Rearranging Equation A. 13 and solving for C, we get 

2 


C = 


‘WvV + 

«> t +p v ) 


WVt 1 

<V p v>‘ 


where 


2 - 3A (P ^~ P v ) 


'o 2w (P £ +P v ) 


2^ 

A(f V P v } 


(A. 13) 


(A. 14) 


(A. 15) 


Mote that C q is the wave speed of a graveling ripple wave in a 

quiet fluid where u = u, = 0. 

v £ 

From Equation A. 14 we can see that the flow will be unstable and 
waves will be formed if 


P 0 P 


< VV > c 2 


*■ v (p t +p v> ° 


or rearranging 


ssasaRW 


u 


P„P 


v£ 


*- V * P £ +P v 


> C 


(A. 16) 


The magnitude of C given by Equation A. 14 varies with A as is shown 
in Figure 5.3. The value of A for which C is a minimum may be designated 
A ^ and corresponds to the point of onset for ripple waves. This value 
of A , can be determined by differentiating equation A. 15 with respect 
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to X and setting it equal to zero. 
Solving for X^ we get 


X . 
rip 


2tt 


y (, v p v ,9 c 


The corresponding minimum value of C is 

o 


o rip 


2°g c | p t -p v ) 

( V°v )2 


1/2 


(A. 17) 


(A. 18) 


For a flowing vapor to cause waves on a liquid surface its speed must be 
greater than or equal to the minimum ripple wave speed, that is, > C q ^ 

We are interested in waves which are more unstable than the weak 

ripple waves, we are therefore interested in waves which have a smaller 

wave length than ripple waves, that is, we are interested in waves with 

X $ ^ r ip' Foj " Sma H values of X the gravity term in equation A. 15 for 

C may be neglected, and then solving for X from Equation A. 16 we get 
o 


X > -2*2— (1 + -s ) 
p vV 1 


(a. 19) 


which is identical to Equation 4.3 given by Tippets [17] ! * 
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